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ABSTRACT 


Tworelated methods for imptoving the dependent f 

*c “Sods in a NASTRAN execution are described. Dan. are patented that 
verify the methods and the NASTRAN DMAP ALTERS. 


1.0 


INTRODUCTION 


a waqtr an user is faced with two major challenges when solving a dynamic eigenvalue problem, 
structure's elastic properties than are required to define ns inertial properties. 

reducing the problem size to one that is more ^ a J^lefo ^ociatetfwith the redistribution of mass 
solution time. 

The theoretical development of the improvement methods is pmvidedin Senior '2- Secti w 
SZ" Conclusions and 

recommendations are provided in Section 5. 


2.0 


THE IMPROVEMENT METHOD 


We begin by deriving the Guyan reduction scheme. 

The dynamic eigenvalue problem is given by the equation 

([K] - X[M]){<|>} = 0 


where 

K = the structural stiffness matrix 
M = the structural mass matrix 
X = the system eigenvalue 
4 > = the eigenvector or modal displacements. 
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Wc can partition Equation 1 into independent DOF, designated in NASTRAN as the analysis set, or A-set, 
and dependent DOF, designated as the OMIT set, or O-set After performing this operation Equation 1 
becomes 
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where die subscript "a" denotes A-set DOF and the subscript "o" denotes O-set DOF. 


Looking at the lower partition of Equation 2 we can say 


( 2 ) 


K L^ + K«oV^LV*M.A-0 (3) 

The Guyan reduction method (Reference 1) makes the assumption that the inertial forces on the O-set 
displacements are much less important than the elastic forces transmitted by the A-set displacements. A 
constraint equation for Guyan reduction can be derived by ignoring the mass terms in Equation 3. The 
resulting constraint equation is given by 


»o=Go1>. (4) 

where 

G o - - i CkL, (5) 

This relationship constitutes a Ritz transformation of the eigenvalue problem. The transformation written 
in terms of the full displacement set is 


«t>) = 



a =[G]{<U 



Using this Ritz transformation, the reduced mass and stiffness matrices become 


( 6 ) 


[M J = [G] t [M][G] (7) 

and 

TKJ = [G] t [K][G] (8) 

The mass of the system is redistributed based upon the elastic connections between the O-set DOF and the 
A-set DOF as shown in Equation 7. 

The reduced mass and stiffness matrices shown in Equations 7 and 8, are then used to compute the 
eigenvalues and the A-set displacements of the reduced system. Once the A-set displacements have been 
computed, the Guyan reduction transformation of Equation 4 is used to recover the O-set displacements. 
This back transformation ignores the inertial terms of the O-set displacements. 
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( 9 ) 


An improved back transformation for * 0 can be foand using Equation 3 (see Reference 2). For mode /. 
this back transformation is given by 

(♦„). - -Poo - W'^ 1 . + 

^^^Sti^^Wroako foe O-se. 4lacemen. recover more efficient. 

Recasting Equation 3 for all the computed modes, we get 

KLO. + K^-MLo.i-M^i-O 00) 

where X is a square matrix with the system eigenvalues along the diagonal. Solving for the ♦» 
displacements that are not multiplied by X, we get 

G 0 <|> t + l + K 0 ^M oo <|) o X = 4> 0 


(ID 


From Equation 1 1 we can see that a closed form solution for fc> does not exist. It is possible, however, to 
use Equation 1 1 to obtain an improved approximation to $ 0 - 

A first approximation to ♦„ can be determined by using the O-set displacements recovered by Equation 4, 

( 12 ) 


or 


Substituting these O-set displacements into Equation 1 1 yields 

G 0 4> a + K>L* # l + h = C 


(13) 


a® are the corrected O-set displacements. These corrected displacements can be substituted back 

into Equation 13 for C and a better approximation, $ 0 , can be computed. This process can be repeated 
SSla~°at the (/ + 1) iteration are the same as the displacement at I 
"super” improved displacements will be identical to those computed using Bluauon 9. and can be 
SSa foe computational penalty associated wifo inverting an O-set by O-se. sued mamx for 

each mode. 

To summarize, force methods for revering foe O-se. fosp^m^rer poforofotg Guyan reduction 

and the reduced eigenvalue analysis have been presented . These three methods are. 

1) Standard Guyan reduction recovery using Equation 4, henceforth designated as Guyan 

2) ta^So-ret displacement recovery using Equations 12 and 13, henceforth designated as 

3) 'SsrivetCSproved O-se. displacements using Equation 13, henceforth designated as 
"super improved" displacements. 
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The reader will note that the A-set displacements are identical for all three methods described above. It is 
assumed that the eigenvalues and the A-sct displacements computed during the eigenvalue analysis are 
"accurate". In other words, the accuracy (or inaccuracy) of the Guyan reduction itself is not in question. 

Thus far, we have discussed improvements only in die O-set displacements. More importantly, any 
quantity computed using these O-set displacements, such as dement forces or dement stresses, will also 
be improved by methods 2 and 3. 

The theory and methodology for improving the O-set displacements has been provided. The following 
section describes the implementation of the improved displacement recoveries in NASTRAN. 

3.0 IMPLEMENTATION IN NASTRAN 

With the methodology in hand, the implementation in NASTRAN becomes an exercise in defining the data 
blocks and the NASTRAN DMAP modules required to perform the desired operations. The DMAP 
ALTER sequences used to recover the improved displacements are provided in Figure 1. The first ALTER 
places the UPARTN module following the SMP 1 module while the second ALTER places the DMAP 
modules used to recover the improved displacements after the SDR1 module. The user controls the 
recovery method with the parameters defined in the DMAP ALTERS. The allowable parameter values and 
the resulting action taken are provided in Table 1 . Note that if no A-set is defined, the O-set recovery 
section is skipped. 


5 

$ DMAP Alter to obtain required matrices for irtprovement . Place after the SMP2 Module. 
ALTER 11 $ where 11 - DMAP statement mmber of Module SMP 2 

UPARTN USET,MFF/,MAOT, ,M00/*F*/*A*/*0* $ 

$ 

5 DMAP Alter to perform O-set displacement improvement. Place after the SDR1 Module. 
ALTER $ where jj - DMAP statement number of Module SDR1 
COND 5KIPIM, OMIT S 

$ 

5 This PARAM defines whether Guyan recovery or inprovement 
S recovery is to be performed (NOIMP < 0, Guyan recovery) 

PARAM / /*NOP*/NOIMP - -1 $ 

COND SKI PIM, NOIMP 5 

S 

$ This PARAM defines what recovery improvement will be performed 
S If NREPT - 0, improve once, NREPT > 0, iterate NREPT times 
PARAM //*NOP*/NREPT -10 $ 

5 

$ MATGEN creates a square matrix from the LAMA table 
MATGEN LAMA/HLAMA/3/2 5 

MPYAD GO, PHIA, /PHIO/O/l/O/ S 

FBS LOO, ,MAOT/Cl/l/l $ 

SMPYAD Cl , PHIA, MLAMA, , , PHIO/A/3///1 5 

FBS LOO, ,MOO/B/l/l 5 

LABEL IMPRV S 

SMPYAD B, PHIO, MLAMA, , , /C/3/ //I 5 

ADD A, C/PHIO/ (1.0, 0.0) / (1.0, 0.0) 5 

REPT IMPRV, NREPT $ 

UMERCE USET, PHIA, PHI0/PHIF/*FV*AV*0* S 

UMERGE USET, PHIF,/PHIN/*N*/*F*/*5* 5 

MPYAD GM, PHIN, /PHIM/0/1/0/ $ 

UMERGE USET, PHIN, PHIM/PHIG/ # G*/*N*/*M* S 

LABEL SKI PM S 


Figure 1. O-set Displacement Improvement DMAP ALTERS 


11 





I Table 1 . DMAP Parameter Settings 1 

KS 

NOIMP 

NREPT 

1 No A-set 

N/A 

N/A 

MM 

-1 

N/A 

wpmm 

0 

0 


0 

# repetitions 


Once the O-set displacements have been recovered, the rest of P' ^ 
This allows the user to define all data recoveries using the familiar NASTRAN (^Control owe 

T^sr'iv-e mems , element forces, element stresses, or any other user requested datawdibe 
JrimJdand handed in the normal fashion. No special provisions are required to view the unproved data. 

4 0 METHOD VERIFICATION 

Two sample problems were created to verify the method and the DMAP describ^xn Soction 3. The first 
sample problem consists of a simple four story building. This problem was used to vwfy the^ 
methodology and the DMAP ALTERS shown in Figure 1. The second problem consists of a 3600 DOF 
substroctuild model. Element forces for this model were recovered from at trmsiem response andysis 
using the three O-set displacement recovery methods and compared to the benchi^de^nt f^es 
obtained when no Guyan reduction was performed. These sample problems verify the improvement 
methods and the DMAP ALTERS. 




m j= 2.0 

k,= 400.0 


m 2 = 2.0 

k 2 = 800.0 


m 3 = 2.0 

kj= 1200.0 


m 4 = 2.0 




k 4 = 1600.0 


K = 


r 400 -400 0 0 

-400 1200 -800 0 
0 -800 2000 -1200 
0 0 -1200 2800 


M = 


r 2 o o o n 

0 2 0 0 

0 0 2 0 

L 0 0 0 2 

r 1/3 2/3 1 

J ° L 0 3/7 J 


f x = 1.278 Hz 
A-set DOF = Uj and u 3 


Figure 2. Simplified Four Story Building 
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The four story building used for sample problem 1 is shown in Figure 2. This problem was selected 
because it is easily represented with NASTRAN elements and may be solved using the NASTRAN 
program. It may also easily be solved by hand so that the data produced by the DMAP ALTERS can be 
verified. Data were recovered for the first mode only. 

Table 2 presents the O-set displacements for the three methods as well as the unreduced benchmark 
displacements. The data in Table 2 were recovered from NASTRAN using the DMAP ALTERS described 
in Section 3. The reader can easily verify that the Guyan results are identical to those recovered by hand 
using Equation 4, the improved results are identical to those recovered by hand using Equation 11, and the 
super improved data are identical to those recovered by hand using Equation 9. These data verify the 
DMAP sequence described in Section 3. 


Table 2. Displacement Comparison 

Disp. 

u i 

Guyan 

Improved 

Super # 
Improved 

Benchmark 
(No A-set) 

1.0000 

1.0000 

1.0000 

1.0000 

u 2 

0.6015 

0.6681 

0.6764 

0.6775 

u 3 

0.4023 

0.4023 

0.4023 

0.4069 | 

U 4 

0.1724 

0.1806 

0.1810 

0.1828 

MAC 

0.99730 

0.99995 

0.99998 

N/A 

* These c 

lata were recovered using 10 iterations 


The Modal Assurance Criterion (MAC) defined in Reference 4 is used to measure the accuracy of the 
eigenvectors provided in Table 2. MAC values will vary between zero, indicating no correlation between 
modes, to unity, indicating perfect correlation between modes. Based on the MAC values, it is clear that 
both improvement methods produce better O-set displacements than the standard Guyan recovery method 
produces alone. 

The advantage of using the improved O-set recovery methods is clearer when element data, e.g. element 
forces or stresses, are compared. The modal spring forces for all three O-set displacement recovery 
methods are compared to the benchmark data in Table 3. From this it is clear that the improved 
displacements produce spring forces that are vasdy superior to those of Guyan reduction. 

Based on this simple problem, the displacements can be dramatically improved by using the methods 
described in Section 2. The next sample problem will show this more clearly. 

The second sample problem uses the 3600 DOF Spacelab Pallet model shown in Figure 3. A transient 
response analysis was performed with this model in two configurations, an unreduced configuration and a 
Guyan reduced configuration. Transient element forces of all die bar elements were recovered using four 
distinct PHASE3 executions, i.e. no A-set, Guyan, improved, and super improved. 
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Figure 3. NASTRAN Model for Sample Problem 2 
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The max imum absolute values for all of the bar forces for the three recovery methods were compared to 
the benchmark case. These comparisons are shown in Table 4. The data are arranged according to a 
percentage difference range. For each of the recovery methods, the percentage of the forces falling within 
this range as well as the maximum difference between the benchmark data and the data produced by the 
three recovery methods within this delta percentage range are provided. 

For example, in the range between two and five percent, 8.53 percent of the forces from the Guyan 
recovery method fell within this range with the maximum difference between the Guyan recovered data 
and the benchmark data being 397. For the improved recovery method, only 0.10 percent of the forces 
fell into this range with a maximum difference between the benchmark and the improved data being 5. The 
percentage of items falling in this range for the super improved method was 0.09, with a m a ximum delta 
of 7. 


Table 4. Bar Element Force Comparisons 
for Sample Problem 2 

A% 

Range 

Guyan 

Improved 


Percentage 
in Range 

Maximum 

|A| 

Percentage 
in Range 

Maximum 

|A| 

Percentage 
in Range 

Maximum 

|A| 

0-2 

89.05 

1045 

99.90 

102 

99.76 

114 

2-5 

8.53 

397 

0.10 

5 

0.09 

7 

5-10 

1.48 

48 

0.00 

0 

0.03 

4 

10-25 

0.60 

82 

0.00 

0 

0.00 

0 

25-50 

0.03 

0 

0.00 

0 

0.00 

0 

>50 

0.32 

2281 

0.00 

0 

0.13 

36 


* These data were recovered using 10 iterations 


As was the case for the simplified model used for sample problem 1, the improved recovery methods 
produce data that are superior to those data computed using Guyan reduction. The data appear to be the 
most accurate for the simple improvement method. This is especially true when the computer CPU time 
required to produce the data is considered. The improved displacement recoveries required 30 percent 
more CPU time than the Guyan recovery, while the super improved displacement recoveries required 150 
percent more CPU time than the Guyan recovery. 

Because of the simplicity of this model, however, it would be premature to use these data to cast the super 
improved method aside without first considering more complex models with equally complex loading. 

5.0 CONCLUSIONS AND RECOMMENDATIONS 

Two methods for improving the O-set displacements were provided. It was demonstrated that both 
improvement methods produce O-set displacements that are vastly superior to those produced using the 
standard Guyan recovery alone. In addition, the NASTRAN DMAP ALTERS required to perform these 
operations were presented along with the supporting data used to verify them. It remains only to 
determine whether the additional accuracy that may be obtainable through the iterative procedure of Method 
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3 is justified by the extra 
guaranteed by the initial 


computational effort After all, a significant degree of approximation is already 
use of Guyan reduction to determine the A-set displacements. 


Because this study did not provide enough information to determine which of the two improved recovery 

methods was bes^ suited for the problems encountered in most engineenng^pbcauonsjit is 

that additional studies be performed to compare improved displacements from a set of models with varying 

complexity to the benchnS unreduced data. These additional cases can ^ 

which improvement method is better in terms of accuracy and computational effiaenc^. RnaUy, « would 

beof nxatintercst to compare the results of a multi-mode transient response analysis before and after 

modalimprovement to assess its practical significance in terms of the end result 


6.0 
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